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Abstract  Third  harmonic  generation  of  a  Gaussian  short  pulse  laser  in  a  tunnel  ionizing 
plasma  is  investigated.  A  Gaussian  short  pulse  laser  propagating  through  a  tunnel  ionizing 
plasma  generates  third  harmonic  wave.  Inhomogeneity  of  the  electric  field  along  the  wavefront 
of  the  fundamental  laser  pulse  causes  more  ionization  along  the  axis  of  propagation  while 
less  ionization  off  axis,  leading  to  strong  density  gradient  with  its  maximum  on  the  axis 
of  propagation.  The  medium  acts  like  a  diverging  lens  and  pulse  defocuses  strongly.  The 
normalized  third  harmonic  amplitude  varies  periodically  with  the  distance  with  successive 
maxima  acquiring  lower  value.  The  self-defocusing  of  the  fundamental  laser  pulse  decays  the 
intensity  of  the  third  harmonic  pulse. 
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1.  INTRODUCTION 


Harmonic  generation  of  electromagnetic  radiation  in  plasma  has  been  a  subject 
of  extensive  study  for  quite  sometime  [1-5].  The  highly  nonlinear  nature  of  the 
interaction  of  the  laser  pulse  with  plasmas  implies  that  harmonic  light  should  be  a 
significant  feature  of  such  interaction.  Third  harmonic  generation  is  a  very  useful 
technique  that  can  convert  output  if  infrared  lasers  to  shorter  wavelength  in  the 
visible  and  near  ultraviolet.  Siders  et  al.[ 4]  have  demonstrated  the  generation  of 
a  tunable  blue-shifted  third  harmonic  generated  at  an  ultra  fast  ionization  front 
with  efficiencies  and  spectral  characteristic  consistent  with  Vanin  et  al.  [5], 
Pressure  tuning  of  the  power  spectrum  centroid  over  the  range  10”5  <  Aco  /  CO  <  10  1 
,  corresponding  to  the  pressure  range  (10  <  p  <  500 Torr).  In  a  nonlinear  medium, 
an  intense  beam  of  finite  diameter,  which  converge  to  a  focus  and  then  diverges, 
may  exhibit  nonlinear  dynamics  that  significantly  affect  propagation  dynamics. 
The  greater  the  intensity,  the  bigger  the  difference  between  incoming  and  outgoing 
beams.  Third  harmonic  generation  with  such  an  input  beam  produces  third  harmonic 
in  the  far  field  when  the  phase  mismatch  between  the  fundamental  and  its  third 
harmonic  is  zero  or  negative  (and  small)  [6].  They  have  also  demonstrated  this 
quantitatively  for  several  cases. 
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Rax  and  Fisch  [7]  have  investigated  relativistic  third  harmonic  generation  in  a  plasma. 
The  phase  velocity  mismatch  results  in  a  third  harmonic  amplitude  saturation  and 
oscillation.  In  order  to  overcome  this  saturation,  they  described  a  phase  matching 
scheme  based  on  a  resonant  density  modulation.  Fedotov  et  al.  [8]  have  observed 
third  harmonic  generation  f  unamplified  30-fs  pulses  of  Cr:  forsterite  laser  in  a  holey 
fiber. 

They  have  demonstrated  that  high  degree  of  light  localization  in  the  core  of  a  holey 
fiber  enhances  third  order  nonlinear  processes,  allowing  the  third  harmonic  of  non 
amplified  30  fs  0.2  nJ  Cr:  forsterite  laser  pulses  appeared  spectrally  broadened  at  the 
output  of  the  fiber.  Lamprecht  et  al.  [9]  have  shown  that  measuring  the  THG-ACF  is  a 
suitable  method  for  evaluating  the  time  characteristics  of  fs  driven  plasmon  oscillations 
in  metal  nanoparticles.  Resonant  and  off-resonant  excitation  of  particle  plasmons  can 
be  clearly  distinguished  by  the  ACF  measurements. 

Ganeev  et  al.  [10]  have  studied  harmonic  generation  from  solid  surfaces  irradiated  by 
27  ps  Nd:glass  laser  pulses  in  the  intensity  range  1013— 1015  W/cm2.  harmonic  emission 
up  to  fourth  order  is  observed  in  the  specular  reflection  direction  with  conversion 
efficiencies  of  2  x  10'8,  10-10  and  5  x  1 0'12  for  second,  third,  and  fourth  harmonic, 
respectively,  for  a  p-polarized  pump  beam  at  1 015  W/cm2.  Sodha  et  al.  [11]  have  studied 
third  harmonic  generation  in  a  collisional  plasma  by  a  Gaussian  electromagnetic 
beam. 

It  is  seen  that  the  self-focusing  of  the  main  beam  enhances  the  power  of  harmonic 
output  by  about  three  orders  of  magnitude  in  a  typical  case  in  the  paraxial-ray 
approximation. 

In  this  paper,  we  study  the  third  harmonic  generation  of  a  Gaussian  short  pulse  laser 
in  a  tunnel  ionizing  plasma.  An  intense  short  pulse  laser  with  a  Gaussian  radial  profile, 
propagating  through  a  tunnel  ionizing  plasma.  The  plasma  has  radial  as  well  as  axial 
inhomogeneity.  The  ionization  starts  on  the  laser  axis  and  creates  a  strong  electron 
density  gradient  with  its  maximum  on  the  axis  of  propagation.  The  medium  acts  like 
a  diverging  lens  and  pulse  defocuses  strongly.  Laser  imparts  an  oscillatory  velocity 
to  electrons  and  exerts  a  ponderomotive  force  on  them,  produces  third  harmonic 
nonlinear  current  density,  driving  third  harmonic  generation.  The  normalized  third 
harmonic  amplitude  varies  periodically  with  the  distance  with  successive  maxima 
acquiring  lower  value.  The  self-defocusing  of  the  fundamental  laser  pulse  decays  the 
intensity  of  the  third  harmonic  pulse. 

In  section  II,  we  obtain  an  expression  for  the  third  harmonic  transverse  current  density, 
deduce  an  equation  for  the  third  harmonic  amplitude  and  solve  this  equation.  A 
discussion  of  results  is  given  in  section  III. 
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2.  SELF-DEFOCUSING  OF  THE  MAIN  PUFSE  AND  NONFINEAR  THIRD 
HARMONIC  CURRENT  DENSITY 

Consider  the  propagation  of  intense  Gaussian  laser  short  pulse  in  a  tunnel  ionizing 
plasma  along  z  direction. 


El  =  M1(z,f)exp[-i(co1t-£1z)] 


(1) 
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fl (z) 
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where  rQ  is  the  spot  size  of  the  main  laser  pulse  at  z  =  0,  and  /  is  the  beam  width 
parameter  of  the  main  laser  pulse.  The  laser  pump  and  the  third  harmonic 
electromagnetic  wave  obey  the  linear  dispersion  relation,  k2  ~  (co2  /  c2)(l-  CO2  /  CO2) 
,  where  C0p  =  (4nn0e2  /  m)m  is  the  plasma  frequency,  c  is  the  velocity  of  light  in 
vacuum,  n0  is  the  electron  density  and  -e  and  m  are  the  charge  and  mass  of  an  electron 
respectively. 


The  evolution  of  plasma  frequency  up  varies  with  time  as12 
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where  0)2(z)  —  4nne2/m  ,  n,„  is  the  initial  density  of  neutral  atoms, 

7  =  (tt/2)1/2(/0  / h)(\E\/ Ea)2  exp(-EA /\e\),  Ea  =  (4 / 3)V2m/03/2  ,  where  I0 
is  the  ionization  potential  and  h  =  lull  is  the  plank’s  constant.  Following  the 
Liu  and  Tripathi13,  for  paraxial  ray  approximation,  we  introduce  dimension  less 
quantities  =  r|  =  ymt'  ,  where  y00=(I0 /h)(n/2)m(l/ g2)  ,  express 

Aj  =  A|()  exp(/LS'J  and  expand  relevant  quantities,  in  the  paraxial  ray  approximation, 
as  S  =  S0  +  S2r2 ,  and  co2  =  co20  +  co22  r2 ,  where  S0  and  S2  are  the  functions  of  z 
and  of  o  =  of  (t, z,r  =  0).  The  equations  for  evolution  of  plasma  frequency  C0/;  and 
beam  width  parameter  fx  are 
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(5) 


where  G  =  exp[(g  -  gl)  /  ff]  ,  gt  =  gf\  ,  g  =  EA/Em  .  In  Eqs.  (3)  and  (4),  G 


characterizes  the  tunnel  ionization  rate  and  fall  off  rapidly  when  beam  width  parameter 
rises.  At  the  front  of  the  laser  pulse  ( 77  =  0),  the  beam  would  diverge  only  due  to 
diffraction. 

The  laser  imparts  an  oscillatory  velocity  to  electrons. 


eE, 
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v, 
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Vj  and  Bi  beat  to  exert  a  ponderomotive  force  =  —(e  /  2c)v,  X  B1  ,  on  the  electrons 
at  20l>[  ,  giving  an  oscillatory  velocity  V2  at  2 lj  , 


It  also  produces  density  perturbation  nl  at  2(0j ,  in  compliance  with  the  equation  of 
continuity.  We  get, 


Density  perturbation  nx  couples  with  Vj  to  produce  nonlinear  current  density  at  third 
harmonic  J,KL  =  --nlevl  . 
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The  nonlinear  current  produces  third  harmonic  wave  with  electric  field  E,  ■  This  field 
also  produces  a  self-consistent  current  J  1 , 


(10) 
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The  wave  equation  for  the  second  harmonic  field  E.t  can  be  deduced  from  Maxwell’s 
equations, 


4k  dJ3  1  d2E, 
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Let  the  complementary  solution  of  this  equation  is 

En  =  A3  exp{-7l3S3}exp[-7'(3co17-&3z)], 


(13) 


where  S3  is  a  function  of  r  and  z  .  Using  Eqs.(9)  and  (10),  and  equating  real  and 
imaginary  parts 
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In  the  paraxial  ray  approximation,  r2  «  r2f3  ,  weexpand  S3  as  S3  =  4>(z)  +  [Lr:  /  2 
,  where  (33  represents  the  radius  of  curvatures  of  the  wave  front  of  the  third  harmonic 
wave.  For  an  initially  Gaussian  beam,  we  may  write 


A2  = 
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(16) 


Substituting  the  values  of  S-,  and  A3  in  Eqs.  (1 1)  and  (12)  and  equating  the  coefficients 
of  r2  on  both  sides,  we  obtain 
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Now  we  suppose  the  solution  (particular  integral)  of  Eq.(12)  is 

E31  =  A'  expf-fXScOj?  -  3&jZ)], 


(19) 


where  A^  =  A'0(z)\|/3,  ip3  =  exp[— r2  /  r02/,2]exp[— 3  ik3S3] .  Using  Eq.(12)  and  (19), 
we  obtain 
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Figure  1:  Variation  of  beam  width  parameter  f  as  a  function 
of  normalized  distance  of  propagation  £  for  different  values  of 
T|  ,  the  other  parameters  are:  eAw  /  m COjC  =  0.02  ,  CO pmrQ  /  c  ~  20 
and  kxc  /  C0j  =  0.44  . 
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Multiply  above  equation  by  \|/3  r  dr  and  integrate  with  respect  to  r  ,  we  get 
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Where  A3^  =  A30  /  A10  ,  we  have  solved  coupled  Eqs.  (3),  (4),  (5)  and  (21)  numerically 
by  applying  the  initial  conditions;  at  77  =  0  ,  co20  =  0  ,  0)(;2  =0,  for  all  at  ;  =  0, 
f  =  1  ,  df  /  dq  =  0  for  all  77.  The  parameters  are  EA  /  Em  =  3  corresponding  to 
I0  =21eV  and  laser  intensity  /  =  6  X 1017  W  /  cm2.  We  have  plotted  in  Fig.  1,  /’  Vs  c 
or  different  values  of  77,  the  other  parameters  are  eA10  /  mcOjC  =  0.02 ,  COpmr0  /  c  ~  20  and 
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kt  c  /  C0j  =  0.44  .  One  may  note  that  at  higher  value  of  T|,  self-defocusing  of  main  laser 
pulse  increases.  In  Figs.  (2)  and  (3),  we  have  plotted  normalized  amplitude  of  the  third 
harmonic  wave  A()  with  normalized  propagation  distance  £,  for  different  values  of  r).  It 


5 

Figure  2:  Variation  of  normalized  third  harmonic  amplitude  as  a  function  of  normalized 
distance  of  propagation  ^  for  f|  .  The  other  parameters  are:  ^AJ0  /  =  0.02 

,  CO  rn  /  C  ~  20  and  k.c  /  (0.  =  0.44  . 
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Figure  3:  Variation  of  normalized  third  harmonic  amplitude  as  a  function  of 
normalized  distance  of  propagation  ^  for  t|  =  1 .  The  other  parameters  are: 
eAn  /  m(0,c  =  0.02 .  CO  r.  /  c  ~  20  and  k.c  /  CO,  =  0.44. 
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is  noted  that  as  the  main  laser  pulse  defocuses,  the  normalized  third  harmonic  amplitude 
decreases  with  for  stronger  defocusing  the  third  harmonic  amplitude  saturates. 

3.  DISCUSSION 
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Third  harmonic  generation  of  a  Gaussian  short  pulse  laser  in  a  tunnel  ionizing  plasma 
is  studied.  An  intense  Gaussian  laser  beam  propagating  through  a  tunnel  ionizing 
plasma  gets  self-defocused  due  to  ionization  induced  refraction.  A  considerable  decay 
in  the  intensity  of  third  harmonic  wave  is  obtained  on  account  of  self-defocusing  of 
the  fundamental  wave.  The  third  harmonic  amplitude  decrease  initially  and  at  the  later 
time  it  saturates.  The  normalized  third  harmonic  amplitude  varies  periodically  with 
the  distance  with  successive  maxima  acquiring  lower  value  and  saturate  for  stronger 
defocusing.  The  self-defocusing  of  the  fundamental  laser  pulse  decays  the  intensity  of 
the  third  harmonic  pulse. 
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